Topological term in the non-linear a model of 
the SO (5) spin chains 



Dung-Hai Lee 

Department of Physics, University of California at Berkeley, Berkeley, CA 94720, 

USA; 

Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, 

CA 94720, USA 

Guang-Ming Zhang 

Department of Physics, Tsinghua University, Beijing 100084, China 

Tao Xiang 

Institute of Physics and Institute of Theoretical Physics, Chinese Academy of 
Sciences, Beijing 100190, China. 



Abstract 

We show that there is a topological (Berry phase) term in the non-linear a model 
description of the SO(5) spin chain. It distinguishes the linear and projective repre- 
sentations of the SO(5) symmetry group, in exact analogy to the well-known #-term 
of the SO(3) spin chain. The presence of the topological term is due to the fact that 
7F2 ( gO(3)xSO(2) ) = ^' ^ e discuss the implication of our results on the spectra of 
the SO (5) spin chain, and connect it with a recent solvable SO (5) spin model which 
exhibits valence bond solid ground state and edge degeneracy. 
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1 Introduction 



The effects of topological terms on the dynamics of Goldstone modes and the quan- 
tum number of solitons and instantons in non-linear a (NLcr) models have a long 
history, and continue to attract strong interests from the physics community[T]. 
For example, in one spatial dimension (ID), quantum SO(3) spin chains have 
fundamentally different low energy properties, depending on whether the site rep- 
resentation is linear (spin integer) or projective (spin half-odd integer). For nearest 
neighbor isotropic exchange interaction, the former cases always have excitation 
gaps while the latter are gapless - the well known Haldane conjecture [2]. Aside from 
the usual stiffness terms, the (1+1)D NLcr models for these spin chains contain a 
topological (Berry phase) term (also known as the 9 term) [2]. When the space-time 
configuration of the Neel order parameter wraps the target space[3] {S 2 ) n times, 
the Berry phase factor is +1 for integer spin chains, while it is (— l) n for half-odd 
integer spin chains [2]. 

Using an algebraic approach, Chen, Gu and Wen[l] recently generalized an idea 
of Ref.[5] and argued that, in one spatial dimension, a gapped ground state which 
is invariant under translation and the global symmetry operation (we refer to this 
type of state as "totally symmetric" in the following) is obtainable when the site 
representation of the global symmetry group (which can be discrete) is linear. 
When the site representation is projective, a totally symmetric ground state must 
be gapless. In a projective representation D, the matrix product D(gi)D(g2), where 
gi t 2 are group elements, can differ from D(gi.g 2 ) by a phase factor e iely9im \ For 
the SO (3) group, an integer spin forms the linear representation while a half- 
odd integer spin forms the projective representation. Hence the spectral difference 
between the translational invariant integer and half-odd integer SO (3) spin chains 
constitutes a special example of the results of Ref.jl]. Thus there are two ways to 
view the difference between integer and half-odd integer SO (3) spin chain: one is 
geometrical (the Berry's phase) [2] and the other is algebraic[l]. 

SO (5) is a rank-2 classical Lie group. It also has linear and projective representa- 
tions. For example, in the vector representation, the generators of the Lie algebra 
are given by [5] 



where a, b, = 1, ..5, and z, j = 1..5. Two consecutive it rotations generated by, e.g., 
L 12 give 



{L a b)jk — —i&a,j&b,k + i&a,k$b,j 



(1) 



U 12 (n)U 12 (n) = I 5x5 . 



(2) 
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The spinor representation, on the other hand, is given by[BJ 



L ab = z[r a ,r b ]/4, (3) 

where T a ^ are the 4x4 gamma matrices (e.g., Ti = — a y <g> cr x , T 2 = —<y y <S> <J y , 
T 3 = —a y <S> cr z , r 4 = a x (g> ^2x2, r 5 = a z <E> 12x2)- In this case it is simple to check 
that two consecutive n rotations generated by L i2 yield 

MtOMtt) = -hxA- (4) 

Thus the vector representation is linear, while the spinor representation is pro- 
jective. According to Ref. [1], a spin chain of the former type can have a totally 
symmetric ground state with a gapped spectrum, while a spin chain of the latter 
type has to be gapless if it is totally symmetric. In the following, we will seek for 
the geometric (Berry' phase) difference between the two cases. 

Another motivation for us to study the Berry's phase of the SO (5) spin chain 
is a recent exactly solvable ID SO (5) spin model (in the vector representation) 
proposed by Tu et. al.[7]. The ground state is a translational invariant matrix 
product state, i.e., a valence bond solid state [8], and the excitation spectrum has 
a gap. Moreover, when the chain is subjected to open boundary condition, there 
are edge states. These properties are reminiscent of those of integer SO (3) spin 
chains [9]. 

With the advance of cold atom physics, the SO (5) spin chain might not be a 
purely academic model anymore. An SO (5) symmetric spin chain can in principle 
be realized experimentally when the hyperfine spin-3/2 cold fermions on an ID 
optical lattice form the Mott-insulating state [TU]. At quarter filling (one fermion 
per site), the effective spin chain is in the spinor representation, while for half- 
filling (two fermions per site) it is in the SO (5) vector representation. Therefore 
the idea presented here might one day be tested experimentally. 



2 Model formulation 



Let us start by considering the following SO (5) invariant Hamiltonian 



K a<b 



\a<b 



(5) 
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where L ab s are the S0(5) generators and Ji j2 > 0. When L ab are given by Eq.((T]), 
the ground state is translational invariant and the spectrum has a gap in the 
parameter range 1/9 < J2/J1 < 1/3 [7] . Naively, one would not expect the NLcr 
model action of this model to contain a topological term. This is because in contrast 
to SO(3) spin chain where the space-time dimension (1 + 1) matches the dimension 
of the target space of the order parameter (S 2 ), for the SO(5) spin chain the target 
space dimension is much larger than the space-time dimension. 

To understand the structure of the target space for the SO (5) spin chain, we need 
to know how the presence of SO (5) "magnetic" moment breaks the global SO (5) 
symmetry. For that purpose, it is sufficient to consider the following mean-field 
theory where non-linear terms in L ab are decoupled into linear ones with order 
parameter given by (L ab ) = {—l) l m ab 

H MF = (-2J X + 2J 2 A 2 ) £ (-l)*m ab 4b + £ (j.A 2 - 3J 2 A 4 ) , (6) 

i,a<b i 

where A 2 = J2 a <b m a b- The question at hand is, for a fixed total magnitude of 
rriab (i.e. fixed J2 a <b m ab)i what is the most energetically favorable ratio between 
different components of m ab . This can be answered by diagonalizing a single-site 
Hamiltonian 

Hx = - E m ab L ab . (7) 

a<b 

and see what ratio gives the lowest ground state energy. (Of course we need to 
remember the sign of m ab change from site to site.) 

In the following, we study the two irreducible representations given by Eqs.([T]) 
and ([3]). First we consider the vector representation, Eq.flTJ. It is straightforward 
to show the energy spectrum of H\ is E = (— A 1; — A 2 , 0, A 2 , Ai) where A 12 = 

sjA± t/W+B-C and 
^ = E™y2 = A 2 /2, 

a<b 

B = 2 ( m acm a drn bc m bd -m ab m ad m bc m cd + m ab m ac m bd m cd ), 

a<b<c<d 

C = E E ml b ml d (l-5 bc )(l-5 bd ). (8) 

a<b a<c<d 

The single site ground state energy reaches the minimum when B = C, where 
the energy spectrum of Hi is E = {—A, 0, 0, 0, A}. Now we ask how is the SO(5) 
symmetry broken under such a condition. Let us use Hi as one of the two Cartan 
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Fig. 1. (color on-line) The root (upper row) and weight (lower row) patterns of the vector 
(a) and spinor (b) representation of SO(5). The x and y coordinates of each dot corre- 
sponds to the eigenvalue of H\ and H 2 . The red arrows indicate how the rasing/lowering 
operators in the SO(5) Lie algebra change the eigenvalues of H\ and f^-The central dots 
in the root patterns are doubly degenerate. 

generators [6]. For the other Cartan generators H 2 , we choose it to be a different 
linear combination of so that Tr(HiH 2 ) = 0. The root and weight diagrams 
are shown in Fig. (JT^l) . Here the x and y coordinates of the dots are the eigen- 
values of Hi and H 2 , respectively. The arrows in the root diagrams indicate how 
the rasing/lowering operators[6] in the SO(5) Lie algebra change the eigenval- 
ues of Hi 2 . H 2 and the rasing/lowering operators generate an SO(3) subgroup 
which commutes with Hi. As the result, the SO(5) symmetry is broken down to 
SO(3)xSO(2), where the SO(2) is generated by H x itself. 



Second we consider the spinor representation, Eq.Q. It is straightforward to 
show that the energy spectrum of Hi is E = (— Ai, — A 2 , A 2 , Ai) with Ai )2 = 

J A ± VC -B/y/2. In this case the single site ground state energy reaches the min- 
imum when A = y/C — B, where the energy spectrum of Hi is E = {—A, 0, 0, A}. 
The root and weight patterns are shown in Fig. ([T])(b). Again the SO (5) symmetry 
is broken down to SO(3)xSO(2). 

After fixing the ratio of m a &, we assume that the low energy fluctuations corre- 
spond to smooth space-time dependent SO (5) rotations of such an order parameter 
pattern. The NLcr model precisely describes such smooth fluctuations. The order 
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parameter lives on the manifold so(S)xSO(2) > wn i cn nas ^ as its second homotopy 
group [TT]. Thus the corresponding NLa model may contain a topological term, 
which can lead to a spectral difference between the vector and spinor representa- 
tions. 



3 The single-site Berry's phase 



To study the possible topological term, we begin by analyzing the Berry's phase 
of a single SO (5) spin described by the following time-dependent Hamiltonian 

H x {t) = -Y,m ah {t)L ab , (9) 

a<b 

where m a b(t) satisfy the constraints: (a) J2 a <b m ab = 1> an d (b) Hi(t) possesses 
SO(3)xSO(2) symmetry. Both constraints can be satisfied by starting with a ref- 
erence Hamiltonian satisfying (a) and (b) and perform time- dependent SO (5) 
conjugation, i.e., 

H 1 (t) = U\t)H lfl U(t). (10) 

As usual, the Berry's phase is given by the loop integral of the Berry connection[TJ. 
We can use Stoke's theorem to convert this loop integral to an areal integral over 
a disk with the loop as the boundary. The advantage of doing so is the Berry 
curvature rather than the Berry connection appears in the latter integral. This 
makes the integral involving only local quantities and gauge invariant: 

S R = f 1 du I dt e" u TrF lIv . 



~2 Jq u J u fc 

In the above F^ v = (d^A v — d u A fl ) and = — i{f)|c? M |fi). Here \Q(t,u)) is the 
ground state of 

H^t, u) = U\t, u)H lfi U(t, u). (12) 

In Eq. (1T2"]) . U(t,u) is the extension of the U(t) in Eq. (ITU]) to the disk. Because 
7r 1 (SO(5)/SO(3)xSO(2))=0, we can always construct the extension so that U(u = 
l,t) = U(t) and U{u = 0, t) = Uq, where Uq is a certain reference SO(5) element. 

Using the first order perturbation theory for wave functions, it is simple to show 
that 

= (n^H\k) { k\d„H\n) 

k (E -E k ) 2 1 ^ 
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Here k labels the excited states. Since all Hamiltonians described by Eq. (I12p are 
unitary conjugate of one another, they have the same eigenspectrum. Under that 
condition we have 

d,H = Y,{E k - E )(\d,k)(k\ + \k)(d,k\). (14) 

k 



In writing down the above equation we have made a shift of the zero of energy 
so that Eq — > 0. Substitute Eq. (Till) into Eq. ([TBI and use the fact that (Q\k) = 
(k\tt) = we find 

TrFp f = -iTi(Q[d lt ,Q,d v Q}), (15) 
where Q is the ground state projection operator 

Q(t, u) = \Q) = U(t, u)PU ] {t, u), (16) 

where P = |0)(0| is the ground state projector operator of H 10 . Substituting 
Eq. (HH into Eq. §IT$ we obtain 

S B = fduj dtTr(Q[d u Q, d t Q\). (17) 

Eq. (TT71) actually applies for any target space. For example, in the case of SO (3) /SO (2), 
we have 

U(t, u)=\ Zl Z2 \ and P = 1 ° 
-z 2 zj loO 



where z±2(t, u) satisfy 

{zx,z 2 ) - a- 

Substitute the above two expressions into Eq. (1161) and Eq. (11 7p we obtain 



h(t, u). 



S B = - [ du [ dt(n- d u h x d t h) , (18) 
2 jo J 

which is the well known expression for the Berry's phase of a spin-l/2[2]. 

Because the dimension of our target space is eight, there are more than one disks 
having the loop in question as the boundary. Therefore, we need to ask whether 
Eq. fTT7|) yields the same answer for the Berry phase when different extensions of 
U(t) — > U (t, u) are used. The difference in the Berry phase using two different 
disks as extension can be calculated by integrating the Berry curvature over a 
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closed two dimensional surface formed by joining the two disks at their common 
boundary. The resulting closed surface has the topology of a 2-sphere. Because 
the second homotopy group of our target space is Z, all 2-spheres in the target 
space are topologically the multiple of a basic sphere. Hence all we need to check 
is whether the Berry curvature integral is integer multiple of 2n when t and u in 
Eq. (I17p parameterize the basic 2-sphere |12j. In the following, we perform such a 
calculation. 

For the vector representation, we choose H^q = L 12 , and pick U(t, u) so that 

U\t, u)H 1)0 U(t, u) = w(t, u) ■ L, (19) 
where L = (L 13 , L 23 , L 12 ), 

Wi(t, u) = (sin(7ru) cos — — , sin(7ra) sin — — , cos(7ra)), 

and j3 is the period of the imaginary time. For the spinor representation, we take 
#i,o = (L12 - L 34 )/y/2, and L = (L 13 + L 24 , L u - L 23 , L 12 - L 3A ) / 

Using Eq. (|T7|) . or equivalently Eqs.flTTJ and (|T3|) . we find that 

C vector _ a q spinor _ cy (OC\\ 

B, basic sphere ) ^B, basic sphere " V I 

Thus the condition for the uniqueness of the Berry phase is satisfied. As we shall 
see, the difference between the vector and spinor Berry phase in Eq. (|2"U|) serves to 
distinguish the vector and spinor SO (5) chains. 



4 The lattice Berry's phase, gapful versus gapless, and the edge states 



We now extend the above single-site Berry phase analysis to the one dimensional 
lattice. If the order parameter is perfectly "antiferromagnetic" at all time, the 
ground state of H\(t) are \Q)(t) on the even sublattice and its conjugate \&{t)) = 
R\fl(t))* on the odd sublattice, where R is the operator that satisfies RL* ab R~ l = 
—L ab . For the vector representation, R is given by R i5 = — i? 2 4 = -R33 = —R42 = 
-R51 = —1 and Rij = otherwise. For the spinor representation, R = —il 2 x2 <E> Qy 
Using the above result, it is straightforward to show that (Q\L a t,\Q) = — (f2|L a &|f2), 
and RU*R~ 1 = U for all U G SO(5). This, plus the invariance of the trace under 
matrix transposition, allows one to show that e^ u TiQd^Qd u Q = —e^ ur YiQd^Qd u Q. 
As the result, the Berry's phases associated with neighboring sites tend to cancel 
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each other. Let r label the center of mass position of site i and i + 1 for i = odd, 
the total lattice Berry's phase is equal to 

^ = E E f du I dtTr(Q r+e/2 [d u Q r+e/2 , d t Q r+e/2 }), (21) 

r e =±l 

where Q is a smooth function of spacial coordinates. Under such a condition, 
Eq. f l2T|) has a continuum limit 

^ ot = \j dx J dtTrQ[d t Q, d x Q], (22) 

the factor of 1/2 arises from the density of odd lattice sites. Under open boundary 
condition, Eq. (122]) becomes 

S T = \j dxdtTiQ[d t ,Qd x Q] 

+ \[du J dt{Tr[Q[d u ,Qd t Q]] R -Tr[Q[d u Q,d t Q}} L }, (23) 

where the subscript "R" and "L" labels the right and the left ends. This topological 
term together with the stiffness term from the energetics, constitute the NLcr model 
for the SO (5) spin chain. 

Eqs. (120]) and (122]) have important implications. The fact that the mapping from 
the space-time to the target space is classified by integer homotopy classes implies 
that the space-time order parameter configurations can be grouped into different 
topological sectors labeled by an integer topological invariant. This is similar to the 
SO(3) NLcr model where the topological invariant, the Pontryagin index [T3]. is the 
number of times which the order parameter configurations cover the target space 
S 2 . In our case there is an analogous integer topological index, which we will refer 
to as the Pontryagin index as well. Eq. ([2"U]) implies that for the vector SO (5) spin 
chain the Berry phase associated with the order parameter configuration having 
different Pontryagin indices are all the same because exp(i47r/2 x integer) = +1. 
Given the facts that (i) the topological term has no effect (hence the NLcr model 
has only the stiffness terms), and (ii) the target space dimension is high, it is easy 
to believe that the vector SO (5) spin chain should have a quantum disordered, 
i.e., translational invariant gapped, phase. For the spinor SO (5) chain, however, 
the order parameter configurations with even Pontryagin index have the Berry 
phase exp(z27r/2 x even integer) = +1, while those with odd Pontryagin index 
have Berry phase exp(z27r/2 x odd integer) = — 1. This result is exactly analogous 
to the Berry's phase in the spin-1/2 representation of the SO(3) antiferromagnetic 
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Heisenberg chain. In view of the result of Ref.jl], we conclude that the above 
non-trivial Berry's phase also implies the lack of an energy gap as long as the 
translation symmetry is unbroken. 

Now we comment on the edge state of the SO (5) "valence bond solid state" in 
Ref. [7]. By tuning the ratio of J\ and J 2 in Eq.©, Tu et al were able to show 
that a short-range entangled, translational invariant matrix product state is the 
exact ground state. In addition, under the open boundary condition the ground 
state wavefunction becomes 4 x 4 = 16 fold degeneracy. According to Eq. (1231) . the 
boundary of a vector spin chain should exhibit the following Berry phase 

^jT 1 duj dtTrQ[d u ,Qd t Q}. (24) 

When Q(t,u) is a unit Pontryagin index order parameter configuration, the value 
of Eq. (l24l) is |*4?t = 2tc. This is consistent with the spinor Berry phase in Eq. ( 120]) . 
Therefore, the edge state of the vector SO (5) spin chain carries the spinor repre- 
sentation. Because the latter is 4- dimensional, each end of the chain independently 
yields a 4-fold degeneracy of the ground state, resulting in a total 4 x 4 = 16 fold 
degenerate ground state for the open chain. 

Before closing a technical remark is in order. The readers might wonder what if 
the ratio between different components of m a b fluctuates away from the the optimal 
value. When that happens the single site spectrum will become {— Ai, — A 2 , 0, A 2 , Ai} 
and {— Ai, — A 2 , A 2 , Ai} for the vector and spinor representation, respectively. In 
this case, the SO(5) symmetry is broken down to SO(2)xSO(2). The second ho- 
motopy group of so ^^so{2) is ^ © ^ rather than Z. In other words, the image 
of the space-time in the target space is topologically the multiple of two basic 
spheres. As A 2 — > 0, one of these spheres shrinks to a point. We have checked that 
so long as A 2 is small, i.e., when the ratio between m a b does not deviate from the 
optimal value too drastically, the Berry phase is only sensitive to the Pontryagin 
index of the dominant (large) sphere. Hence all the results discussed earlier remain 
unchanged. 



5 Conclusion 

We have studied the Berry's phase of the antiferromagnetic SO (5) spin chain, and 
shown the existence of a topological term in the non-linear a model description 
of the system. The quantum phase factor associated with this topological term 
differentiates the vector (linear) and spinor (projective) representations. We argue 
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that this leads to the spectral difference as long as the translation and SO (5) 
symmetry is unbroken. More specifically, the vector spin chain can have a totally 
symmetric ground state while having an energy gap. The spinor chain, on the other 
hand, must be gapless if there is no symmetry breaking. Under the open boundary 
condition, we find the boundary Berry's phase of the vector spin chain is consistent 
with the derived edge degeneracy of an exactly solvable model. The present result 
can be straightforwardly generalized to other irreducible representations, leading 
to two classes of SO (5) spin chain: in one class the site representation is linear, and 
in the other the site representation is projective. While the first class can have a 
totally symmetric ground state and maintain a gapped spectrum, the second class 
must have gapless spectrum if there is no symmetry breaking. 

Finally, we would like to make contact with relevant works in the literature. 
Eq. ([T7|) appeared in Ref.fH] in discussing a super-symmetric version of the Hub- 
bard model. Of course, the target space is entirely different. In Ref.[15] Affleck 
discussed the topological term of an SU(N) spin chains. In that case the target 
space is U(N)/U(m) x U(N — m), which also has Z as its second homotopy group. 
In Ref . [TE] . Read and Sachdev have discussed the effects of Haldane's (SU(2) or 
SO(3)) Berry phase term in two space dimensions [M]. Subsequently, they gener- 
alized this to the SU(N) group in Ref.[T7]. The essential difference between the 
SU(N) and the SO (5) spin chains is their target spaces. Different target spaces 
have distinct topology. As far as we know, this is the first time one points out 
that the target space of the SO (5) spin chain also has second homotopy group Z, 
and constructs the Berry's phase term for the NLu model description. Of course, 
analogous to Refs. [TIIfTT] one can ask what is the effect of this topological term in 
two space dimensions. We leave this question for future studies. 
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